Design, biological evaluation, and molecular modelling insights of cupressic acid derivatives as promising anti-inflammatory agents

Abstract The major labdanes in the oleogum resin of Araucaria heterophylla (Salisb.) Franco, 13-epi-cupressic acid (1) and acetyl-13-epi-cupressic acid (2) were used to prepare seven new (3–9), along with one known (10) derivatives. RAW264.7 cells were used to evaluate the anti-inflammatory activity of the derivatives (1–10) via measuring the level of COX-2 expression and IL-6. Pre-treated RAW264.7 cells with 1–10 (except for derivative 7) at 25 µM for 24h exhibited downregulation of COX-2 expression in response to LPS stimulation. Moreover, pre-treatment with compounds 1, 2, or 3 significantly attenuated the LPS-stimulated IL-6 level in RAW264.7 cells (p < 0.05). A docking study was conducted against phospholipase A2 (PLA2), a crucial enzyme in initiating the inflammatory cascade. The significant structural features of compounds (1–10) as PLA2 inhibitors included the carbonyl group at C-4 (free or substituted) and the hydrophobic diterpenoid skeleton. This study suggested 13-epi-cupressic acid as a scaffold for new anti-inflammatory agents.


Introduction
Since ancient times, plant-derived products have been used to reduce inflammation or inflammation-associated conditions 1 . Araucaria trees are coniferous trees grown as ornamental plants, which upon injury, exudate an oleogum resin that protects against invading pathogens 2 . Few studies reported the biological activities of the oleogum resin of Araucaria heterophylla; mainly as an antiulcerogenic, cytotoxic, and antibacterial agent 3 . The oleogum resins, containing mainly labdane diterpenes, have been employed topically to treat bruises and dermal inflammation 4 . Diterpenes are natural products known to exert potent anti-inflammatory activities through various mechanisms 4 . Labdane-diterpenes were reported to exhibit anti-inflammatory activities in several animal and cell-based bioassays 3,4 . In addition, the analgesic and antiinflammatory effects of semisynthetic terpene derivatives were previously described 4 . Cupressic acid derivatives are the major labdane diterpenes of the oleogum resin of Araucaria heterophylla. Previous studies by our research group reported the antiprotozoal activity of resin-isolated cupressic acids and the anti-inflammatory activities of their microbially transformed metabolites 2,5 .
Physiologically, inflammatory reactions are normal immune responses against potential or already existing tissue injury. However, pathologically it may lead to, or append, multiple complications in many diseases, such as cardiovascular, pulmonary, and rheumatic diseases; in addition to malignancies, and skin disorders 6 . Phospholipase A2 (PLA2) plays a crucial role in initiating the inflammatory cascade in mammalian cells 7 . Inflammatory stimuli disrupt the cell membrane, and PLA2 starts its action on the phospholipid bilayer and releases arachidonic acid (AA), a polyunsaturated fatty acid 8,9 . At this stage, overexpression of cyclooxygenases (COXs) and lipoxygenases (LOXs) arises due to high concentrations of their AA substrate. Successively, this leads to the formation of various inflammatory mediators, including cytokine interleukins, ILs (such as IL-6 and tumour necrosis factoralpha, TNF-a), prostaglandins, PGs (such as PGE2 and PGD2), thromboxanes, TXAs (such as TXA2, TXA4, and TXB4), leukotrienes, LTs (such as LTA4, LTB4, and LTC4), and lipoxins, LXs (such as LXA4 and LXB4) [8][9][10] . Therefore, targeting PLA2 may lead to the discovery of new PLA2 inhibitors, which can mitigate inflammation by down-inhibiting these inflammatory reactions (Figure 1).
Little has been done on the semi-synthesis and biological evaluation of labdane terpenes obtained from Araucaria heterophylla. Cupressic acid diterpenes (Figure 2), possess a carboxylic acid moiety and a hydroxyl functional group, both of which are amenable to chemical modification and functionalization. In this work, we implemented the approach of rational, yet simple, chemical derivatization to semi-synthesize more potent and/or selective anti-inflammatory lead compounds. Additionally, a docking study has been conducted to explore the possible inhibition of PLA2 by the investigated compounds.

General
The 1 H and 13 C Nuclear Magnetic Resonance (NMR) spectra were measured in CD 3 OD, with tetramethyl-silane (TMS) as the internal standard, on BRUKER Ascend TM (400 MHz) spectrometer (Bremen, Germany) operating at (400 MHz) for 1 H and 100 MHz for 13 C NMR. Chemical shifts (d) were reported relative to TMS, as part per millions (ppm); coupling constants (J) in (Hz), while multiplicities of signals were expressed as: (s ¼ singlet, brs ¼ broad singlet, d¼ doublet, dd ¼ doublet of doublet, m ¼ multiplet). Processing of NMR data was done using MestReNova, v. 12.0.1-20560 (Mestrelab Research S.L., Santiago de Compostela, Spain). IR spectra were reported using a Thermo Scientific Nicolet TM iSTM10 FT-IR spectrophotometer (WI, USA). For reaction monitoring, thin layer chromatography (TLC), with precoated Silica gel aluminium plates 60 GF 254 (Merck, Nagel, Germany) was used. Any other used solvents or reagents were of commercial grade.

Plant materials
The resin exudates from the trunk of Araucaria heterophylla (Salisb.) Franco trees (Araucariaceae) were collected from the Gardens of Mansoura University campus, Mansoura, Egypt, in February 2021. The plant identity was confirmed by staff members of the Horticulture Department, Faculty of Agriculture, Mansoura University, Egypt.

Molecular docking
The in vitro assays as well as in silico target prediction studies implied that the mode of actions of prepared derivatives (1-10) might be mimicking that of corticosteroids, which encouraged us to run molecular docking studies to support the findings.

Protein preparation
The PDB protein crystal was processed by protonation, bond corrections, and fixing the potential using the MOE software. Also, the bound H 2 O molecules were removed. The active binding pocket was specified using the site-finder tool in the MOE software and it was selected by isolating atoms and backbone "dummies" at alpha centres on the basis of the 3D coordinates of the inhibitor's nearby residues.

Preparation of the studied ligands
The structures of the studied analogs were sketched using the ChemDraw Professional 15.0.0.0.106 # 1998-2015 (PerkinElmer Informatics Inc., MA, USA). This was followed by protonation, assigning partial charges by the current forcefield tool, and minimising energies by the general energy minimisation tool.

Docking experiment
Molecular docking was performed by the Triangle Matcher for Placement and Rigid Receptor for Refinement. London DG and GBVI/WSA DG were selected as the scoring functions. The binding free energy of the studied analogs was calculated as kcal/mol. The highest number of negative scores obtained was chosen as the highest ligand-enzyme interaction score for the studied structures.
Physicochemical and ADME (absorption, distribution, metabolism, and excretion) predictions The physicochemical and ADME properties, in addition to druglikeness, and toxicity were in silico investigated for the synthesised Figure 3. Summarised schemes (i-v) of the general preparation procedures of cupressic acid derivatives (1-10). (i) 1.5 equiv. of acetic anhydride and 2 equiv. of pyridine at room temperature for 24h; (ii) in DCM, reflux for 3h with 1.5 equiv. thionyl chloride (SOCl 2 ) or for 6h with an excess of SOCl 2 ; (ii-a) in DCM, reflux with 2 equiv. triethylamine (TEA) and 1.5 equiv. ethanolamine; (ii-b) in DCM, reflux with 2 equiv. 2-phenyl ethylamine at 45 C for 48h; (iii) reflux in dry acetone, excess potassium carbonate, and a catalytic amount of DMAP; (iv) Stirring with excess selenium dioxide solution (SeO 2 ) in anhydrous dioxane. analogs using the SwissADME online tool (Swiss Institute of Bioinformatics; http://www.sib.swiss) 12 .
Evaluation of the anti-inflammatory activity of the tested compounds RAW264.7 cells were treated either with 5 mg/mL of LPS (Escherichia coli (O111:B4), L2630, Sigma-Aldrich) or pre-treated with test compounds at 25 mg/mL for 3h then treated with 5 mg/mL of LPS for 24h. After that, the culture media was collected for determination of the IL-6 level by an ELISA kit (FineTest TM , Q-M0121-B, Wuhan, China), and the cells were collected for total RNA extraction by TRIzol reagent (Invitrogen TM , Cat. No. 15596026, Carlsbad, CA, USA). Cells pre-treated with 1 mM of dexamethasone (as an anti-inflammatory drug), then treated with LPS, were used as the positive control 13 .
RNA concentration and purity were evaluated by a NanoDrop (Thermo Scientific TM , Waltham, MA, USA). RNA (4 mg) was reverse transcribed to cDNA by using Top script TM RT DryMix (dN18/dN6 plus, Cat. No. RT220, Enzynomics, Daejeon, Republic of Korea) kit. The level of COX-2 and b-actin was detected by SYBR Kit (SensiFAST TM , Meridian Bioscience, Cincinnati, OH, USA) by using qPCR (StepOne TM Real-Time PCR, Applied Biosystems TM , Waltham, MA, USA). The expression level was first normalised to the level of b-actin and then calculated as 2 À(DDCt) . Primers (Eurofins Scientific#, Alachua, FL, USA) are listed in Table 1.

Statistical analysis
The statistical analysis was performed by comparison of means obtained from triplicate treatments with the One-Way Analysis of variance "ANOVA" method on GraphPad Prism Version 8.0 (GraphPad Software, Inc., San Diego, CA). Differences between groups were significant at "p values" of <0.05.

Results and discussion
Preparation of the proposed derivatives (1-10) The starting material (13-epi-cupressic acid, Torulosic acid, 1) was isolated from the oleogum resin of Araucaria heterophylla (Salisb.) Franco and subjected to acetylation reaction to yield acetyl-13epi-cupressic acid, 2 (Scheme i, Figure 3). A reactive acyl chloride derivative was produced from the protected acetyl compound (2) and reacted with a nucleophile (ethanolamine) to produce the amide derivative 3 (Scheme ii-a, Figure 3). Interestingly, upon repeating the same reaction under more drastic conditions, oxidative deacetylation took place to form communic acyl chloride which upon reacting with a nucleophile (phenylethylamine) produced the carbamoyl derivative, 4 (Scheme ii-b, Figure 3). The ester derivatives (5-9) were produced through the coupling of the potassium salt of compound 1 with different alkyl halides (Scheme iii, Figure 3). Allylic hydroxylation of the bicyclic diterpene (10) was achieved through the oxidation of compound 1 with SeO 2 (Scheme iv, Figure 3).

Structural identification of the prepared derivatives (1-10)
Spectroscopic (FT-IR, NMR & MS) and elemental analyses were employed to confirm the structures of the prepared compounds, 1-10 ( Figure 3). Free OH/NH groups showed broad absorption peaks at (3600-3200 cm À1 ). Upon inspection of NMR spectra, convergent values -with slight variationswere observed among all analysed compounds (1-10), as expected, due to the common bicyclic system with its exocyclic double bond and carbonyl group in addition to the chain-terminal double bond. The 1 H NMR spectra (Table S2) 13 C NMR spectra (DEPT/APT experiments) displayed extremely downfield quaternary carbon signals at d C 183.2-171.9 (C-19) assigned to the carbonyl carbon next to the bicyclic system (Table S1). The exocyclic double bond showed a downfield quaternary carbon signal at d C 150.9-147.7 (C-8) and up-field secondary carbon at d C 107.7-106.6 (C-17), while the chain terminal double bond showed a downfield primary carbon signal at d C 147.9-141.5 (C-14) and upper-field secondary carbon at d C 113.5-109.8 (C-15). Generally, the upfield tertiary carbon signals, including two at d C 30.4-22.9 (C-16/C-18) and one at d C 13.0-12.4 (C-20), represented the common 3 methyl groups. Other eleven upfield aliphatic carbon signals (Table S1) indicated the rest of the typical backbone of labdane-diterpene in all derivatives. Compounds 2 and 3 showed extra methyl signals corresponding to the methyl groups of the introduced acyl moiety, in each case. This includes the resonances at d H 1.96 (-COCH 3 , in 2 and 3). Also, the 13 C NMR spectra of 2 and 3 showed an extra downfield quaternary signal at d C 169.8-169.3 (-COCH 3 ) and an upfield methyl signal at d C 22.1-20.9 (-COCH 3 ) corresponding to the acetyl group on the side chain-hydroxyl group at C-13.

Characterisation of compounds 3 and 4
The FT-IR max (neat) spectra of amide derivatives (3 and 4) showed absorption bands for (NH) and (C¼O) functional groups in the regions t (3400-3000 cm À1 ) and (1750-1710 cm À1 ), respectively. The 1 H NMR spectrum of compound (3) showed the presence of a methyl singlet at 1.96 of the acetyl group at C-13 (Table S2 and Figure S6). Whereas no acetyl group was detected in the case of 4, it rather showed a downfield triplet at d H 5.36 (H-12) and a relatively upfield multiplet at d H 2.33 (H-11), which indicated the introduction of an additional double bond between C-12:13 Table 1. Primer sequence and their reference in NCBI.

Gene
Primer sequence Reference sequence COX-2 creating a terminal conjugated diene (Table S2 and Figure S13). Due to the newly introduced conjugated double bond, the 13 C NMR and APT spectra of compound 4 (Figures S14, S15) showed two extra downfield carbon signals, including a methine at d C 134.0 (C-12) and quaternary at d C 133.5 (C-13). Further evidence for the loss of the esterified OH group was indicated by the absence of a quaternary signal at d C 83-85 (C-13) compared to the related derivatives (2 and 3 and the presence of upfield methyl carbon signals at d C 11.9 (C-16). Regarding the introduced substitutions at C-19, the 13 C NMR and APT spectra of compound 3 showed two characteristic methylene signals at d C 42.8 (N-CH 2 ) and 61.5 (O-CH 2 ), for C-1 00 and C-2 00 of the ethanolamine moiety, respectively. Similarly, the 13

Characterisation of compounds 10
The FT-IR max (neat) spectra of compound 10 showed prominent peaks for (O-H) and (C ¼ C) groups at t (3446 cm À1 ) and (1692 cm À1 ), respectively. The 1 H NMR spectrum of compound 10 ( Figure S39) showed more downfield split (triplet) signal at d H 4.29 (H-7) corresponding to the proton geminal to the newly introduced hydroxyl group at C-7 (Table S2 and Figure S39). The 13 C NMR and APT spectra of compound 10 ( Figure S40 and S41) showed a more downfield oxymethine signal at d C 74.8 (C-7) corresponding to the newly introduced hydroxyl group located at the allylic carbon to the exocyclic double bond 14 . Thus, the structure of compound 10 was determined as 7a-hydroxy-13-epicupressic acid, which was previously obtained before by microbial transformation of cupressic acid by Fusarium graminearum 14 .

Biological assay
Effect of cytotoxicity of the tested compounds against RAW264.7 macrophage cells RAW264.7 is a known model for the in vitro investigation of antiinflammatory activity 15 . To inspect the safety of the investigated derivatives (1-10) on RAW264.7 macrophage cells, a cytotoxicity study was performed. WST-1 reagent was used to detect the viability of the treated cells. The derivatives were toxic to RAW264.7 at 100 and 50 mg/mL, while more viable cells were observed at 25 mM. Therefore, all the experiments were performed at a safe dose to the cells (i.e., 25 mM).

Modulation of COX-2 expression in RAW264.7 cells
Based on the cytotoxicity assay which showed that the investigated compounds were safe to RAW264.7 cells at 25 mg/mL, the anti-inflammatory activity was evaluated. LPS is an external cell stimulus for immune cells via binding with Toll-Like receptor 4 17 . Stimulation with LPS produces inflammatory cytokines, such as IL-6 and COX-2 enzyme. In this study, LPS treatment led to a change in the macrophage morphology (elongated and spindle shape) while the pre-treatment of macrophage cells either with compounds 1, 2, 3, 4, 5, 6, 8, 9, or 10 or Dexamethasone inhibited the response to LPS and decreased the inflammation (Figure 4). The RAW264.7 cells treated with LPS alone expressed a high level of the inflammatory gene COX-2, while pre-treatment of the cells for 3h with the tested compounds (except compound 7) showed downregulation of COX-2 expression in response to LPS. The results were comparable to the standard anti-inflammatory drug, Dexamethasone (Figure 4). In contrast, compound 7 which is characterised by the presence of a hydrophobic hexyl chain at C-19 of torulosic acid, stimulated COX-2 expression. The obtained results were in full agreement with the previously published reports which demonstrated significant anti-inflammatory activity of different diterpene skeletons on LPS-induced RAW264.7 cells by various mechanisms, including the downregulation of COX-2 gene expression [15][16][17] . The first labdane-type diterpene to be discovered as an anti-inflammatory agent was cis-communic acid isolated from the leaves of Cryptomeria japonica (Taxodiaceae). It was identified as the main biomolecule responsible for the activity of a topical traditionally used Japanese preparation for the treatment of skin eruption, swelling, eczema, and injury. Cis-communic acid was shown to inhibit carrageenan-induced paw edoema in rats and histamine-induced contraction in guinea pig ileum 18,19 .
Determination of the level of IL-6 by ELISA Recognition of LPS by macrophage cells leads to the activation of downstream signalling pathways and transcription factors, such as Nuclear Factor-jB (NF-jB) 20 and Activated Protein-1 (AP-1) 21 . The activated transcription factors induce the expression of different cytokines, including IL-1, IL-6, and TNF-a. The anti-inflammatory activity of the compounds was further confirmed by measuring the level of IL-6. To investigate the effect of derivatives on proinflammatory cytokines, the level of IL-6 was measured in the medium of RAW264.7 cells after treatment either with LPS alone or pre-treatment with the studied compounds for 3h, followed by treatment with LPS for 21h (Table 2 and Figure 5). LPS treatment dramatically enhanced the expression of IL-6. Analysis of the results revealed that pre-treatment with compounds 1, 2, or 3 significantly attenuated the LPS stimulation and reduced the level of IL-6 (p < 0.05). It has been documented that IL-6 plays a major role in the induction of septic shock via releasing prostaglandins  and leukotrienes and activating signal transducer and activator of transcription 3STAT3 22 . The bioactive derivatives 1-3 shared common structural features, including the occurrence of terminal hydroxyl groups in the substituents at C-19 of torulosic acid, in addition to the absence of the hydroxyl group at C-7. This report described for the first time the effect of 13-epi-cupressic acid (Torulosic acid) and its derivatives on the downregulation of inflammatory cytokines expression, such as IL-6. However, several studies found in the literature demonstrated a similar effect exhibited by other diterpene classes [23][24][25] .
Previous reports showed that labdane diterpenoids may exert their anti-inflammatory action through multiple mechanisms mainly found at the levels of gene expression and/or enzyme activity, such as interfering with the NF-jB signalling pathway, NO production, or arachidonic acid (AA) metabolism 26 . Forskolin andrographolide, sclareol, and (þ)-polyalthic acid are examples of labdane diterpenoids reported in the literature as potential anti-inflammatory constituents from several successful traditional remedies 26 . Our results demonstrated that most of the prepared derivatives were able to downregulate COX-2 gene expression and attenuate IL-6 level in LPS-stimulated RAW264.7 cells. In conclusion, our findings suggested 13-epi-cuppressic acid and its derivatives as promising leads for the development of new anti-inflammatory drugs. However, further in vivo experiments are required to study their pharmacokinetic and toxicological properties.

Docking experiments
Docking experiments were performed using MOE software. The docking results of the prepared derivatives (1-10) against the active binding pocket of PLA2 are displayed in Table 3. It demonstrated that the greatest docking score was recorded by derivative 10 (À8.97 kcal/mol), which showed greater binding affinity than that of the co-crystallized inhibitor, BR4 (À8.01 kcal/mol). These are followed by the starting compound 1 (À7.87 kcal/mol), and derivatives 2 (À7.86 kcal/mol). Whereas, the lowest binding affinity was recorded for 7 (À6.34 kcal/mol), followed by 9 (À6.72). The binding affinities are in full agreement with the obtained in vitro assays as compound 7 was not able to downregulate COX-2 gene expression. On the other hand, derivatives 1 and 2 were among the most active derivatives in both Leukocyte accumulation and IL-6 assays.
PLA2 inhibitors are mainly characterised by their ability to form electrostatic interactions with Ca þ2 ion, H-bonding with the catalytic amino acid, His47, and various types of hydrophobic binding interactions with the hydrophobic residues lining the active site of the enzyme 27 . The structure of co-crystallized ligand (BR4) was docked within the investigated set of derivatives for validating the docking results, where it showed remarkable binding affinity (À8.01 kcal/mol) and the same interactions as the crystal structure, including Ca:301 (ion contact), Gly29 (H-bond and arene-H), Gly31 (H-bond), and His47 (H-bond) ( Figure S47). However, Dexamethasone showed the lowest binding affinity (À5.81 kcal/ mol) and no binding interactions ( Figure S3).
Visualisation of the docked structures revealed that all compounds showed interactions with the crucial Ca þ2 metal ion (Ca:301) except for 3 and 7 ( Figure S49 and S50). Also, all compounds did not show interaction with His47 residue. However, most of them showed H-bonding interactions with Gly29, Gly31, and Lys52 11,27 .
Compound 1 showed binding interactions with the Ca:301 (metal contact), residues Gly31 (H-bond), and Lys62 (ion contact) through the COOH group at C-4 ( Figure 6(a)). Similarly, 2 interacted with Ca:301 and Gly31 via the COOH group at C-4 ( Figure  6(b)). The aliphatic ester derivative 5 showed a remarkable capacity to form several interactions mainly via the introduced ester group at C-4, which was able to extend the molecules to interact with Ca:301 (metal complex), Glys31, and Lys62 (H-bonding), in addition to the parent's carbonyl at C-4 that showed H-bonding with Gly29 ( Figure 6(c)). Extended conformation was found to be required for optimal PLA2 inhibitory activity for filling the hydrophobic active site of the enzyme 11 . Remarkably, derivative 10 showed the top score (À8.97 kcal/mol) due to various interactions with the active site of PLA2. These interactions include two Ca:301 metal contacts with both the COOH group at C-4 and the OH group at C-7. Moreover, it exhibited H-bonding interaction with Gly29 through the OH group at C-7 ( Figure 6(d)).
The aromatic ester derivative 9 similarly showed interactions with the parent's carbonyl at C-4 with both Ca:301 (metal contact) and Gly31 (H-bond). But it was able to form an arene-H interaction with Lys52 through the introduced benzyl moiety. The aliphatic ester derivative 6 only showed a metal ion interaction with Ca:301 through its OH group at C-13. Compound 3 also showed one interaction with the hydrophobic amino acid Leu2 through H-bonding with the introduced ethanolamide moiety. However, the less polar hexyl derivative 7 did not show any interactions, which can be further explained by its lowest binding affinity (À6.34 kcal/mol). This result agreed with the obtained in vitro anti-  inflammatory activity for this derivative. The isobutyl ester 8 showed H-bonding with Cyst44 via the OH group at C-13, and the ester group at C-4 formed a metal complex with Ca:301. Finally, derivative 4 showed a metal contact with Ca:301 and H-bonding with Gly29 residue through the parent's carbonyl group.
It could be concluded that the most important structural features required for the activity of 13-epi-cupressic acid derivatives as PLA2 inhibitors include the presence of carbonyl group at C-4 which is crucial for binding to Ca þ2 metal ion as was observed in most compounds (Table 3). Additionally, the hydrophobic diterpenoid skeleton is essential for filling the hydrophobic part of the active site (Figure 7). Also, the presence of aliphatic substitution at COOH (C-19) is more favourable than aromatic ones, although free unsubstituted COOH showed better binding affinities and in vitro anti-inflammatory activity as observed in 1, 2, and 10.

Physicochemical and ADME predictions
The physicochemical, ADME, and drug-likeness properties of the prepared derivatives (1-10) were in silico investigated using the SwissADME online software (Swiss Institute of Bioinformatics; http://www.sib.swiss) 12 . The results (Table S3) showed conformity with the Lipinski's rules. However, one violation was recorded for compounds 4, and 7-9 due to their elevated lipophilicity (MLOGP: Moriguchi octanol-water partition coefficient >4.15), 28 . Likewise, none of the prepared derivatives showed PAINS alerts (Pan-assay interference compounds). Additionally, all derivatives showed high GIT absorption as indicated by the predicted values' range of 0.55-0.85, indicating oral bioavailability. Also, no possible inhibitions were recorded for any of the investigated compounds against CYP450-1A2 indicating that they are safe and less likely to affect the metabolism of drugs 29 . The blood-brain barrier (BBB) represents an obstacle in the delivery of drugs to the brain, derivatives 1-3, 8, and 10 were able to cross BBB but 4-7 and 9 were not 30 . The P-glycoprotein (P-gp) contributes to multidrug resistance to several drugs through the efflux of drugs back into the blood 31 . Compounds 3-5 and 10 may show practical resistance based on the predicted results. Radar map is an interpretation of the average values of six descriptors applied for the quick assessment of drug-likeness properties. These descriptors include the size of the molecule, lipophilicity, polarity, flexibility, saturation, and solubility ( Figure S51). The synthesized derivatives were found to fall inside the pink area of their corresponding radar plots 12 . Interestingly, the bioactive compounds 1, 2, and 10 were located entirely inside this pink area which indicated that they have the properties of drugs.

Conclusions
Diterpenes, including the labdane group, have been reported to possess both in vitro and in vivo anti-inflammatory activities. In this study, we prepared a series of semisynthetic derivatives based on the naturally occurring labdane diterpenoids namely, 13-epicupressic and acetyl-13-epi-cupressic acids obtained from the trunk of Araucaria heterophylla (Salisb.) Franco trees. The chemical reactions used in the preparation of eight derivatives (3-10) from the starting molecules were described, including esterification, amidation, and allylic hydroxylation reactions. The anti-inflammatory activity of compounds was evaluated in RAW264.7 cells via measuring the level of COX-2 expression, and level of IL-6. The results showed that most of the prepared derivative exhibited significant (P < 0.05) anti-inflammatory activities as revealed from downregulation of COX-2 expression, diminished IL-6 level of LPSinduced inflammation in RAW264.7 cells. An in silico docking experiment on phospholipase A2 (PLA2), a vital enzyme in starting the inflammatory cascade, revealed the importance of the carbonyl group at C-4, (free or substituted with aliphatic chains) and the hydrophobic diterpenoid skeleton. This study suggested 13-epi-cupressic acid as a scaffold for the synthesis of new antiinflammatory agents. The current study proposed 13-epi-cupressic acid as a promising natural lead molecules for the discovery and development of new anti-inflammatory drugs.